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We present schemes to prepare two-dimensional cluster states [H. J. Briegel and R. Raussendorf, 
Phys. Rev. Lett. 86, 910 (2001)] with atomic ions confined in a microstructured linear ion trap 
and coupled by an engineered spin-spin interaction. In particular, we show how to prepare a n x 2 
cluster state by creating a linear cluster state and adding third-neighbor entanglement using selective 
recoupling techniques. The scheme is based on the capabilities provided by segmented linear Paul 
traps to confine ions in local potential wells and to separate and transport ions between these wells. 
Furthermore, we consider creating three- and four-qubit cluster states by engineering the coupling 
matrix such that through the periodicity of the time evolution unwanted couplings are canceled. 
All entangling operations are achieved by switching of voltages and currents, and do not require 
interaction with laser light. 

PACS numbers: 03.67.Bg, 03.67.Lx, 37.10.Ty 



I. INTRODUCTION 

Cluster states [T| are the physical resource needed 
for a one-way quantum computer [2J |3] — a scheme for 
measurement-based quantum computation [4]. It has 
been shown that the one-way quantum computer pro- 
vides a universal set of quantum gates. Furthermore, 
cluster states can be used to efficiently simulate quan- 
tum circuits. On the other hand, some quantum circuits 
based on the cluster-state model cannot be interpreted 
as a network of gates such as the bit-reversal gate [5]. 

Not only are cluster states a central ingredient for 
measurement-based quantum computing, they are also 
of interest for investigating questions of fundamental rel- 
evance, for instance, regarding the robustness of entan- 
glement. The presumption that the lifetime of entan- 
glement decreases with the number of constituents of 
an entangled system and that, therefore, entanglement 
does not have to be taken into account when describ- 
ing the properties of mesoscopic and macroscopic sys- 
tems is widespread. This is indeed true, for example, for 
Greenberger-Horne-Zeilinger s6 states for which the de- 
coherence rate increases linearly with the number N of 
qubits. However, such behavior does not necessarily hold 
for other TV-particle systems. In fact, many-body states 
exist where the entanglement between the constituents 
is expected to not decay faster than that determined by 
the decoherence rate of a single qubit j7] . Cluster states 
belong to this class of entangled states. 

The topology of cluster states may exist in different 
dimensions. In [8 it is shown that one-dimensional (ID) 
cluster states can be efficiently simulated by classical 
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computers, thus two-dimensional (2D) cluster states are 
needed for useful quantum computation. The first exper- 
imental realization of cluster states was reported in [5]. 
In that experiment, the entanglement was created using 
controlled collisions between atoms confined in an opti- 
cal lattice. Photonic one-way quantum computers have 
already been used to implement quantum algorithms, 
namely, Grover's algorithm [TO] and the Deutsch algo- 
rithm on a four-photon cluster state |llj . Four- and six- 
qubit cluster states have also been recently experimen- 
tally realized with photons [T2l IT5] . 

Until now, multi-qubit cluster states have not been cre- 
ated with trapped ions. In [TJ] it was proposed to gen- 
erate a linear cluster state of four, five, and six ions in 
a linear Paul trap using a gate similar to the M0lmer- 
S0renson gate [T51 HH] ■ This would be accomplished by 
collective addressing of all ions by means of two laser 
beams tuned to the blue- and red-sideband transition of 
the exploited vibrational mode. This technique is robust 
against heating and does not require the ions to be cooled 
into their motional ground state. However, scaling this 
method to larger numbers of ions becomes complex. 

In this paper we propose schemes to prepare two- 
dimensional cluster states with atomic ions confined in a 
linear (one-dimensional) trap. In contrast to the scheme 
described above, we consider creating cluster states by 
means of spin-spin coupling induced by a magnetic-field 
gradient that creates a state-dependent force acting on 
each qubit. Here, laser light is not required to achieve an 
entangling gate between the ions. Furthermore, we con- 
sider two-dimensional cluster states, and suggest that our 
scheme should be highly scalable. In addition, cooling the 
ion string to its motional ground is not necessary. 

An easily scalable method for preparing cluster states 
is applying a Hamiltonian equivalent to an Ising-type in- 
teraction on qubits, initially prepared in states |+) = 
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The scheme for creating 2D cluster states detailed in 



Sees. II A and II B is underpinned by simulations of elec- 
trostatic potentials in a microstructured ion trap. To be 
(1) concrete, we used for this purpose the parameters of a 
microtrap that is currently being developed. 



The time evolution e ~' i - ff */ fi |+}® results in cluster states, 
if Ja.a't = 7T + 2kn with k € N and restricting the inter- 
action to next neighbors yields two-dimensional cluster 
states. However, in general, experimentally accessible J 
couplings J a . a i are of such a form that preparing cluster 
states (or performing other operations) is not trivial. 

This paper presents methods to engineer spin-spin cou- 
plings suitable for generating two-dimensional cluster 
states, a prerequisite for one-way quantum computing 
with state-of-the-art linear ion traps. The outline of the 
paper is as follows: Sec. |I A| is a review of how one- 
way computing with two-dimensional cluster states can 
be implemented on a smaller number of physical qubits 
by reusing those qubits that have been measured during 
the computational process. 

In Sec. |IB| we summarize the relevant properties of a 
collection of spin-coupled trapped ions and then explic- 
itly show in Sec. [II] how to prepare a n x 2 cluster state 
with a number of operational steps that is linear in n 
using this system. 

Section |II A| deals with schemes for cluster state gen- 
eration, where the local electrostatic trapping potential 
experienced by each ion is individually adjustable. We 
outline two schemes that arc primarily suited for creating 
one-dimensional cluster states, even though they could, 
in principle, be used to generate two-dimensional clus- 
ters. For implementing the first scheme, proposed by Mc 
Hugh and Twamley [17l [18] , the harmonic oscillator fre- 
quencies that characterize the local trapping potentials 
are set to equal values, resulting in uniform nearest neigh- 
bor couplings. Controlled generation of nearest neighbor 
couplings by a choice of appropriate non-uniform trap fre- 
quencies is considered in the second scheme. Then, with 
the help of numerical simulations based on an existing 
microstructured ion trap, we show that these methods 
require smaller trap-electrode structures than currently 
available to attain coupling constants that are useful in 
practice for cluster state generation. 

While the above mentioned schemes are based on 
nearest-neighbor interactions provided by individual po- 
tential wells, Sec. |II A2| deals with the question of 
whether coupling constants can be engineered in such 
a way that they fulfill periodicity conditions imposed on 
the time-evolution operator suitable for generating clus- 
ter states in one time-evolution step. 

In Sec. |IIB| we propose a scalable scheme for creating 
two-dimensional cluster states, which does not rely on a 
static placement of the ions in individual wells. Instead 
it makes use of the possibilities provided by segmented 
traps to adiabatically transport ions, and separate ions 
held in common traps, thus allowing us to circumvent 
problems arising in the previously mentioned schemes. 



Breaking up a n x m cluster state into clusters 
of size n x 2 



In order to make one-way quantum computing possible 
first a two-dimensional n x m cluster state needs to be 
prepared and then single-qubit adaptive measurements 
in qubit-specific bases are performed in order to achieve 
a quantum gate [3]. The basis of a measurement may 
depend on the outcome of previous measurements. 

The question this paper deals with is: "How can a 
two-dimensional cluster-state be efficiently prepared with 
trapped ions?" Generally, preparing a two-dimensional 
cluster-state turns out to be difficult for an arbitrary 
number of qubits. In order to work out an experimentally 
feasible procedure for generating cluster states, we notice 
that - within the scope of the one-way model - a n x m 
cluster state may be broken up into segments of dimen- 
sion n x 2 as shown by Mc Hugh [17] and recapitulated 
for the reader's convenience in what follows. 

The usual way to implement a measurement based 
quantum computer is to prepare the entangled state fol- 
lowed by single-qubit measurements. This procedure 
is equivalent to the following: imagine the qubits ar- 
ranged in a two-dimensional array. Then, the initial 
state containing the input data for the quantum gate 
to be performed is written upon the first column of 
qubits. The second column is prepared in the state 
|+) = ^(|0) + |1)). After entangling all qubits with 
their nearest neighbors via the time evolution generated 
by Hamiltonian (I]) with J a ,a't = ft, & n x 2 cluster state 
is established. Now, measurements on the first column 
are performed in a predetermined basis such that these 
measurements amount to executing a desired quantum 
gate. Due to the entanglement, the quantum information 
generated by these measurements is thus transferred to 
the second column's qubits. At the same time, the en- 
tanglement between the two columns is erased by the 
measurement operation. 

That means, the first column may now be used to sim- 
ulate the third column of a n x m cluster. For this pur- 
pose, its qubits are prepared in the state |+) and entan- 
gled with the second column's qubits. Then, single-qubit 
measurements are performed on the second column so 
that the quantum information is transferred to the first 
column's qubits. Thus, performing this procedure m — 1 
times, one can simulate a n x m cluster via anx2 cluster. 
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B. Spin-coupled trapped ions 



We consider N ions of mass m confined strongly in two 
spatial dimensions (x and y direction) and by a weaker 
potential in the third (z direction). This can be achieved, 
for example, by using a linear clcctrodynamic trap where 
the effective harmonic confinement in the radial direc- 
tion, measured in terms of the secular angular frequency 
v r , is much stronger than in the axial direction charac- 
terized by secular frequency v z [HO E01 HI] . The ions are 
laser cooled such that their radial degrees of freedom are 
frozen out and they form a linear Coulomb crystal. 

Each ion provides two internal levels (e.g., hypcrfinc 
states) serving as qubits described by the Pauli operator 
a z . In addition, a magnetic field B — (Bq + bz)e z is 
applied such that the ions experience a gradient b along 
the z-axis. Then, the Hamiltonian of the system reads 



JY 



N 



H = - ^ U n (zo,n)Vzm + ^ hv n o) n a n (2) 



n=l 



The first term of the Hamiltonian represents the internal 
energies of N qubits where the qubit resonances are given 
by u>n(zo,n)} and 2o,n is the axial equilibrium position of 
the nth ion. The second term expresses the collective 
quantized vibrational motion in the axial direction of the 
ions with eigenfrequency v n of vibrational eigenmode n. 
The last term describes a pairwise spin-spin coupling be- 
tween qubits with the coupling constants. 



2m ^— ' u 2 dz n 

j— i J 

h duj n . dbJ r , 
2 S^' 20 '" lh n 



dz 10 



(3) 



where A is the Hessian of the trap potential and D is the 
unitary transformation matrix that diagonalizes A. The 
eigenvalues of A are given by mv 2 . Therefore, the spin 
chain can be interpreted as an iV-qubit molecule with 
adjustable coupling constants J nm , an ion spin molecule. 
If the ions are confined in a global harmonic potential, 
then J nm oc b 2 jv\. 

In Fig. [T] as a concrete example, the coupling con- 
stants J nm are displayed for eight m Yb + ions in a linear 
trap characterized by v\ = 2ir ■ 200 kHz and a magnetic 
field gradient of 100 T/m. 

So far we have a considered a global axial potential 
confining the ion string. However, the trapping poten- 
tial for each ion may be shaped such that the ions reside 
in individual harmonic potential wells. This is accom- 
plished by dividing the dc electrode of a linear Paul trap 
into segments to which individual voltages may be ap- 
plied that shape the axial potential experienced by the 




FIG. 1: (Color online) Coupling constants for eight Yb + ions 
in a Paul trap with axial trap frequency v\ = 2-7T x 200 kHz 
and a magnetic field gradient of 100 T/m. 



ions. Thus, one or more ions may be held in local poten- 
tial wells, and there are additional handles to change the 
range and strength of the coupling constants J nm |18j . 
Recently, segmented microstructured traps have been in- 
vestigated experimentally and theoretically. Such traps 
provide the capability of storing ions in separate poten- 
tial wells, and of separating and transporting ions into 
different trap regions [24j EH [26j EH EHl ES] . 

The spin-spin coupling mediated by the vibrational 
motion in Eq. ^ arises when the ions are exposed to 
a magnetic field gradient that induces a state-dependent 
force. The scheme for cluster-state preparation proposed 
here can also be applied to the case when the required 
spin-spin coupling is generated by means other than a 
magnetic-field gradient. In |30j it was shown that an op- 
tical state-dependent force may induce a coupling whose 
formal description is identical to what is outlined above. 
Electrons confined in an array of microstructured Pen- 
ning traps and exposed to a spatially varying magnetic 
field also exhibit a similar spin-spin coupling [33J . 



II. PREPARING CLUSTER STATES USING 
SPIN-SPIN INTERACTIONS 

Spin-spin coupling as it appears in Eq. ^ may be used 
to prepare cluster states. This is achieved in two steps 
PQ: first, all qubits are prepared in the state |+). Second, 
the spin-spin coupling according to Eq. Q is switched 
on for a time such that / J nm dt = f + 2fc7r, k G N for all 
qubit pairs (n, m) that are to be entangled. This way of 
preparing cluster states provides, in principle, an efficient 
and scalable method to generate entangled states. 

In actual experiments, the above condition can be ful- 
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FIG. 2: Preparing a 4 x 2 cluster-state: In a first sequence of 
operations, a linear 8-qubit cluster-state is created (indicated 
by solid lines connecting the qubits). The second sequence 
then yields entanglement of third-neighbor qubits (dashed 
lines) resulting in the two-dimensional cluster topology. 



filled by manipulating the coupling constants, by apply- 
ing a pulse sequence to selectively realize specific cou- 
plings, or by combining these two methods. In what 
follows, we will address the issue of creating suitable in- 
teractions for cluster state preparation in detail. 

With a string of trapped ions, if all ions are located in 
the same harmonic potential well, the J couplings vary 
in strength throughout the ion string (see Fig. [TJ. Thus, 
achieving controlled dynamics of the system is possible 
only at high cost, for example by using selective recou- 
pling pulse sequences whose length grows quadratically 
with the qubit number |34j . Furthermore, the vibrational 
motion of the ion string is mediated by the Coulomb re- 
pulsion. So J couplings decrease for non-nearest neigh- 
bors. This turns out to be a serious problem for an effi- 
cient preparation of large two-dimensional cluster-states. 

The coarse procedure for creating nx2 cluster states 
that forms the basis for the schemes presented in the 
remainder of this paper (unless noted otherwise as in 
Sec. Ill A 21) is illustrated in Fig. |2l First, a linear n 



qubit cluster state is prepared. Subsequent entanglement 
of third-neighbor qubits then results in the desired two- 
dimensional cluster. 

Alternatively, one could create two four-qubit linear 
cluster states (e.g., qubits 1 to 4 and qubits 5 to 8). 
These linear graph states could then be converted into 
box cluster states by local operations and relabeling of 
some qubits as described in [35 . Thereafter, entangling 
qubits 4 and 5 and qubits 3 and 6 would result in the 
same graph as depicted in Fig. [2] 



neighbor interactions could be achieved. In that scheme, 
individual harmonic oscillator potentials (characterized 
by trap frequencies typically of order 1 MHz) confine a 
linear array of ions such that they are spatially separated 
by 10 /im. The relatively small distance between neigh- 
boring potential wells is necessary to achieve reasonably 
large coupling constants J nm . 

In this arrangement of trap potentials, nearest- 
neighbor couplings dominate. Second (third)-neighbor 
couplings reach values of w 1/6 (« 1/25) of the nearest- 
neighbor couplings (these specific values result from the 
choice of the two parameters trap frequency and ion 
separation mentioned above). Because of the small 
third-neighbor couplings in this arrangement, the gen- 
eral scheme sketched in Fig. [2] is not well suited for 
cluster state generation. In [T7] it was proposed that 
a 2D cluster state can be created by utilizing nearest- 
and second-neighbor couplings in separate steps. This 
scheme would then require appropriate refocusing pulse 
sequences to eliminate undesired couplings (e.g., third- 
neighbor couplings). 

An alternative scheme using microstructured traps 
with electrode dimensions of order 10 /im or smaller is 
sketched in the following. One may set the potentials of 
these individual traps such that, at a given time, only ions 
i and i + 1 interact. This will be the case, if p[ u\, vl +1 
with j =/= i, i + 1. This choice of the strength of individual 
potential wells ensures strong suppression of non-nearest- 
neighbor couplings (below we give a concrete example). 
Confining the ions in such a trap configuration and ap- 
plying a magnetic field gradient for a suitable time such 
that J J i i+ idt = 7r/2 results in maximal entanglement 
between ions i and i + 1. Thus, a linear cluster state can 
be obtained by subsequently performing this operation 
on ions 1 through N — 1. In order to create a n x 2 clus- 
ter state, third-neighbors need to be entangled, that is, 
ions 1 and 4, 3 and 6, 5 and 8. This may be accomplished 
by first setting v\ through vf to a frequency much lower 
than the remaining frequencies to enable coupling be- 
tween ions 1 and 4. Now, ions 1 and 4 are entangled uti- 
lizing spin-spin coupling and selective refocusing is used 
on ions 1 through 4 to undo the unwanted couplings in 
this quartet of ions, thus realizing the coupling J14. The 
other third-neighbors arc entangled analogously. 

In order to check the feasibility of implementing this 
scheme with currently available ion traps we performed 
numerical calculations of electrostatic potentials achiev- 
able in a typical microstructured trap. 



A. Ions confined in individual potential wells 

The inhomogeneity of spin-spin couplings in the case 
of a single harmonic potential can be treated by locating 
the ions in the individual potential wells created in mi- 
crostructured traps. In |17j it is shown that by placing 
each ion in an individual potential well, uniform nearest- 



1. Realization with microstructured traps 

The schemes for generating cluster states outlined 
above require ion traps with electrode structures at a 
characteristic length scale of around 10 ^m in order to 
achieve coupling constants in the kilohertz regime. Even 
though such structures appear feasible, they imply that 
the distance between ions and a solid-state surface is of 
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FIG. 3: (Color online) Sketch of the geometry of a segmented 
ion trap considered for potential simulations. The electrode 
layers are separated by s = 350 )im, while separation in radial 
direction is g = 250 fim. The thickness of the electrodes is t = 
125 /im, the axial length k = 100 )im, and the axial isolation 
distance is h = 30 /jm (compare [41]). 



the same order of magnitude, which leads to significant 
heating rates of the ions' secular motion [351 |371 I3"B] . 
This, in turn, is likely to impede precise quantum logic 
operations [39] - On the other hand, in recent experi- 
ments strongly reduced heating rates have been observed 
with ion traps in cryogenic environments [37, 38 . Thus, 
by sufficiently cooling ion traps, this difficulty that arises 
with small electrode structures may be overcome. 

For many existing microstructured traps typical axial 
lengths of electrode segments are of order 100 /zm (see 
Fig. [3]), and isolation spacings between electrode seg- 
ments are typically of order 30 /im. Such segmented mi- 
crotraps could serve to create individual harmonic oscilla- 
tor potentials for each ion. However, with such relatively 
large trap structures this would lead to large mutual dis- 
tances between ions and thus to small coupling constants 
J. This will make it difficult to employ the schemes out- 
lined above for efficient cluster state generation as will 
be shown now by means of a concrete example. 

For numerical simulations we used the parameters of 
a microtrap that is currently being developed. This trap 
is a three-layer microstructured segmented trap with two 
trapping regions. The upper and lower layers both carry 
electrodes for applying rf and dc electric fields. The 
middle layer serves as a spacer and contains segments 
of current-carrying coils that generate a spatially vary- 
ing magnetic field [40]. We consider a trapping region 
in our potential simulation with the following geomet- 
ric parameters of the electrodes: the two electrode layers 
are separated by the distance s — 350 /im and the seg- 
mented electrodes are separated in radial direction by 
a gap g — 250 /im. The thickness of the electrodes is 
t = 125 /im, and the axial length of each electrode seg- 
ment amounts to k = 100 /im. For isolation, the electrode 
segments are divided by a gap of h — 30 /im. 

We consider 17 segment electrode pairs for the poten- 
tial simulation. The simulation was created using the 
isim package |41j . which is based on boundary-element 
methods. All coupling constants calculated in this sec- 
tion are done for six 171 Yb + ions. 

This ion trap is structured such that for each ion an in- 
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6 
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Voltage / V 


48 


-8 


48 





37.1 


18.3 


27.4 


18.3 


36.8 


Segment No. 


10 


11 


12 


13 


14 


15 


16 


17 




Voltage / V 





48 





48 





48 


-8 


48 





TABLE I: Example for voltage configuration for the mi- 
crostructured trap described in the text to establish the single 
coupling J23 in a string of six m Yb + ions while strongly sup- 
pressing all other couplings. 



dividual potential well may be applied and thus the trap 
frequencies can be individually set. It is then possible 
to generate a sequence of single nearest-neighbor cou- 
plings with suppressed non-nearest- neighbor couplings. 
As noted above, setting two neighboring oscillators' fre- 
quencies to small values compared to all other's allows 
selective coupling of a single pair of qubits. 

Coupling of, for example, only ions number 2 and 3 can 
be attained by applying the voltage configuration given 
in Table|l]to the trap dc electrodes. These voltages result 
in the simulated potential shown in Fig. [4] A polynomial 
fit up to second order yields trap frequencies for oscilla- 
tors 2 and 3 of 0.35 MHz and 0.27 MHz, whereas all other 
frequencies are between 0.8 and 1.6 MHz. Furthermore, 
the distance between oscillators 2 and 3 is smaller by a 
factor of around 2 compared to the other oscillator dis- 
tances. The nearest-neighbor coupling constants shown 
in Table [H] illustrate that J23 dominates over all other 
couplings by 2 orders of magnitude as desired. But due 
to the large distance between ions of 140 /im, even this 
dominating coupling is very small. So segmented mi- 
crotraps with much smaller axial electrode lengths and 
isolation spacings (of order 10 /im) would be required to 
achieve J couplings in the kilohertz range. 
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FIG. 4: (Color online) Potential simulation for the voltage 
configuration shown in Table [I] The arrows indicate local po- 
tential minima (i.e., the equilibrium position of the respective 
ion). Oscillators 2 and 3 have small frequencies compared to 
the other oscillators so that J23 is the dominating coupling. 
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d,. 


i+i/fim 


Vi /MHz 


Ji^+i/Hz 


1 




320 


1.65 


0.001 


2 




138 


0.35 


0.610 


3 




297 


0.27 


0.006 


4 




266 


1.16 


0.001 


5 




279 


0.83 


0.001 


6 






0.98 





TABLE II: Distances di^+i between the individual poten- 
tial wells i/i + 1 and axial trap frequencies Vi determined 
from polynomial fits up to second order around the potential 
minima shown in Fig. [4] These parameters clearly lead to 
the domination of the coupling J23 by two orders of magni- 
tude over all other couplings. In this way controlled nearest- 
neighbor couplings may be implemented without the need for 
refocusing pulses. To increase J-couplings smaller distances 
between ions, and therefore smaller electrode structures are 
required 



2. Periodicity of the time-evolution operator 



One may also utilize the periodicity relation to realize 
only NN interactions. In this case, the J matrix reads: 

f + 2fc 3 7T 2k 2 Tr 2fci7r 

f + 2fc 4 7T 2k 2 TT 



f + 2k 3 TT 



e - / jdt - 1 4 + 2k37r 1 

2fcl7T 2fc 2 7T I + 2fc 3 7T 

(7) 

with ^1,^2,^3,^4 € N. The time evolution Q describes 
an effective nearest-neighbor interaction, thus generating 
the linear four-qubit cluster state: 



|C 4 ) =exp(z-(tj Zi icr Zi 2 + o- z ,20\j, 3 + Vz,30- z a))\+) 6 
=l. u . ^(|0000) + |0011> + 11100) - |1111)) 



(8) 



Replacing the matrix elements O14 — 2k\i: = O41 in 
matrix [7] by 0i4 = 2fci7r + 7r/4 = O41 would allow for 
generating a 2D cluster state. 



The schemes outlined above for preparing a two- 
dimensional cluster-state are based on the generation of a 
linear cluster-state of ions and subsequent third-neighbor 
couplings. These third-neighbor couplings, while undo- 
ing the unwanted next-neighbor (NN) couplings, could 
be accomplished by selective rccoupling techniques (com- 
pare Sec. II B I. However, simultaneous coupling of all 
qubit pairs would be advantageous. In this section we 
show that tailoring the time-evolution operator, that is, 
imposing a suitable periodicity condition by sculpting the 
J couplings, allows for creation of a cluster state for three 
qubits, which corresponds to a triangle graph, and the 
creation of a linear cluster state for four qubits in one 
time-evolution step. 

The general form of the time-evolution operator is: 



(4) 



i<j 



where Qtj — J Jijdt. In order to obtain cluster states 
by spin-spin coupling, the 0^- need to take on values of 
|,| + 27r, I +47r, whereas 6 values of 0, 27r, 47t, ... 
transform the time-development operator into the iden- 
tity operator. So using periodicity to entangle three 
qubits requires the following J matrix: 

- + 2fc7T — \ 

H= / .1,11= I J +2fc7T 4 f +2fc7T (5) 

\ f+2fc7T ) 

with k G N. Applying this time evolution to qubits pre- 
pared in |+)® results in the three-qubit cluster state 



|C 3 ) = exp(ij(a zA a Z:2 + <t z ,ict z ,3 + c z ,20V3))|+)' 5 
= Lu . *(|000> + |111» . 
Here l.u. denotes equivalence up to local unitaries. 



(6) 




FIG. 5: Locating ions in individual microtraps and superpos- 
ing a harmonic long-range potential serves to create coupling 
constants suitable for periodic entanglement. 



In the following we show how coupling matrices can 
be achieved that fulfill the desired periodicity. As before 
(Sec. ~ 



II A 1 1 , the ions are placed in individual harmonic 
oscillator potential wells with adjustable frequencies. Su- 
perposing a long range harmonic potential affecting all 
ions, one more degree of freedom is available (sketched 
in Fig. [5]). So the problem reduces to finding trap fre- 
quencies resulting in J couplings that fulfill the period- 
icity relation. The coupling constants are functions of 
the ions' equilibrium positions, which can be calculated 
analytically only for two and three ions. Furthermore, 
a change in only one trap parameter, say in one trap 
frequency, affects all couplings simultaneously So due 
to the highly non-linear nature of the coupling constants 
as well as the sensitivity to parameter alterations, find- 
ing trap configurations that are suitable for utilizing the 
periodicity condition becomes increasingly difficult with 
the number of ions involved. Here we present empirically 
found parameters. 

Consider three Yb + ions in individual potential wells 
superposed by a harmonic long-range potential of fre- 
quency uj = 2tt x 100 kHz. Trap frequencies of traps 1 
and 3 are 27r x 277 kHz. The middle individual trap has 
frequency 27rx 100 kHz (see Figs. [5]and[6]). The potential 
wells are separated by a distance of 20 /im, and a mag- 
netic field gradient of 100 T/m is applied along the trap 
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axis. These parameters result in the following J matrix: 

/ 7.85 0.87\ 
J/(100 Hz) = 7.85 7.85 (9) 



These couplings are useful for cluster state preparation 
since the periodicity relation is fulfilled: 



21 



hi 



= 9.02 



2tt + 7r/4 



(10) 



N 




FIG. 6: (Color online) Coupling constants suitable for gen- 
erating a triangle graph state in one time evolution for three 
Yb + ions in individual potential wells superimposed on a har- 
monic long-range potential of frequency lu — 2n ■ 100 kHz. 
Trap frequencies of traps 1 and 3 are 2ir ■ 277 kHz. The mid- 
dle individual trap has frequency 2tt ■ 100 kHz. 

In an analogous fashion, we present parameters use- 
ful for generating a four-qubit NN interaction. The 
four ions are again located in individual traps of fre- 
quency uji = — 2tt x 415 kHz and u 2 = w 3 = 
2ir x 280 kHz. A harmonic long-range potential of fre- 
quency u — 2ir x 239 kHz is then applied. The individual 
traps are separated by 5 /im, and the magnetic-field gra- 
dient is again |^ = 100 T/m. With these parameters, 
we obtain the following J matrix: 



.7/(100 Hz) 



/ 4.33 2.08 1.05N 

4.33 4.36 2.08 

2.08 4.36 4.33 

Vl.05 2.08 4.33 



(11) 



which results in following periodicity relations: 

4 • 2tt + 7r/4 



J32 
J 4 1 
J21 
J41 
J31 
J41 



4.15 



4.12 



2tt 

4 • 2tt + 7r/4 

2k 
2 • 2tt 



(12) 



2tt 



We conclude from the case of three ions that periodic 
entanglement is in principal possible, but as can be seen 
from the case of four ions, finding appropriate trap fre- 
quencies becomes more and more difficult with higher 
number of ions due to the non-linearity of the problem. 
We conjecture that numerical approaches such as genetic 
algorithms could be efficiently applied to find optimal 
parameter configurations. 

In Sec. |II A l| it was demonstrated that suitable voltage 
configurations can be found to implement a desired cou- 
pling matrix. In order to attain reasonably large coupling 
constants small electrode structures (of order 10 pxa) are 
required. For the scheme presented in Sec. II A2[ relying 
on appropriate periodicity of the time-evolution opera- 
tor, again small electrode structures are required to be 
implemented efficiently. 



Including ion transport for generating cluster 
states 



Now we turn to the description of a different scheme 
for generating cluster states, which takes advantage of the 
capabilities of segmented ion traps, even with relatively 
large electrodes. Segmented ion traps may not only serve 
for creating an array of linearly arranged potential wells, 
but also for transporting ions into different trap regions 
as well as for separation of ions held in a common trap 
into two distinct traps [24] , l25 l l26 l l27 l [28 ] [29] . 

Adiabatic ion transport using segmented microtraps 
has been demonstrated by Rowe et al. In this ex- 

periment a 9 Be + ion was transferred between trap loca- 
tions 1.2 mm apart in 50 /is with almost unit efficiency. 
Furthermore, separation of two ions held initially in a 
common trap into distinct traps was demonstrated. This 
was accomplished by using a five electrode configuration. 
Fast non-adiabatic transport of 40 Ca + ions was reported 
in [29 . The experimental results show a success rate of 
99.0(1)% for a transport distance of 2 x 2mm in a round- 
trip time of T = 20 /is. Application of optimal control 
theory is planned to achieve lower excitation of vibra- 
tional motion in the future. 

A n x 2 cluster is prepared using the scheme illustrated 
111 Fig. [2| This generation of cluster states is accom- 
plished in two operational sequences. During the first 
sequence nearest-neighbor couplings are established that 
lead to the preparation of a one-dimensional cluster state. 
The second sequence then establishes couplings between 
third-neighbor qubits and serves to create a 2D cluster 
state of eight ions. These two sequences will be detailed 
in what follows. First, we outline the sequences and state 
the required time evolution. Then, numerical simulations 
of the potentials of a microstructured ion trap, which is 
currently being tested, will serve to illustrate the feasi- 
bility of the proposed scheme. 
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1 1.6 1 0| 2 0| 2| 0| 2 0| 2| 0| 2| 0| 2 0| 2| 0| 1.6 



TABLE III: Voltage configuration for creating 8 harmonic po- 
tential wells such that oscillators 2 to 7 have uniform axial 
frequencies of 200 kHz. Voltage at electrode pairs 1 and 17 is 
lower in order to reduce fringe effects. 



1. First sequence: creating a ID cluster state 

The first sequence itself consists of two steps. During 
the first time step with duration t\ , pairs of ions, namely 
ions number 1 and 2, 3 and 4, 5 and 6, 7 and 8 occupy 
a common trap potential. Switching on the magnetic 
field gradient for a desired time results in N/2 uniform 
nearest-neighbor couplings: 



The duration t\ is chosen such that Jl 1 Jdt 
filled and the time-evolution operator 



is ful- 



U tl = exp(i-(a Zyl a z . 2 



is obtained. 

At the end of the first time interval, the magnetic- 
field gradient is turned off (i.e., the spin-spin is zero), 
and the ions sharing a common trap are separated and 
transported into potential wells such that ions 2 and 3, 
4 and 5, 6 and 7 occupy a common trap. When the 
ion transport is finished, the magnetic-field gradient is 
switched on again during time ti, thus resulting in the 
other half of NN couplings: 

7T 

U t2 = exp(i-(<r z .2<J z ,3 + cr z ,4,cr Zt5 + cr zfi a z j)) (15) 

The time evolution during the first sequence described 
above (i.e., between t = and t = t^) requires N/2 po- 
tential wells of equal axial trap frequency for TV ions. We 
now show that the required axial potential can be realized 
using the segmented microtrap introduced in Sec. |I B| 
Applying a voltage of 1.6 V to the outermost electrodes, 
and an alternating series of 2 V and V, respectively, on 
the remaining dc electrode pairs (Table |III[ ) , results in 
harmonic axial potentials of frequency to = 2n x 200 kHz 
(compare Fig. |7| 

Six of the eight created harmonic oscillators have the 
same frequency, so that the scheme for preparing a two- 
dimensional eight-qubit cluster state can be realized with 
this axial potential. Keeping the pairs of ions 1 and 2 
through 7 and 8 in potential wells 2 through 5 for the 
time ti results in the time evolution according to Eq. ( 14 ) 
when the magnetic-field gradient is switched on. There- 
after, the field gradient is deactivated, and the ions are 
rearranged, such that ions 1 and 8 occupy potential wells 
2 and 6 respectively, and the pairs of ions 2 and 3 through 
6 and 7 are located in wells 3 through 5. Reactivating the 
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FIG. 7: (Color online) Simulated potential for the voltage 
configuration shown in Tab ]III| Six of eight potential wells 
have uniform axial frequencies of 200 kHz as well as equidis- 
tant axial spacings of 260 /jm. 



TABLE IV: Distances between the minima of the potential 
wells in /im for the voltage configuration shown in Tab. |III| 
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260 


260 


260 


260 
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field gradient for a time ti leads to the evolution given 
in Eq. ( 15 1 , thus creating an eight-qubit linear cluster 
state. 



The distances between the minima of the potential 
wells are shown in Table II VI The characteristic value of 
the distances between the traps 2 to 7 is simply given by 
2 • (k + h) = 260 /mi. Deviating values of oscillators 1 and 
8 can be explained by fringe effects. 

Trap frequencies of to = 2-k x 200 kHz and a magnetic- 
field gradient of d z B = 100 T/m result in coupling con- 
stants of « 3 kHz. In Fig. [8] the coupling constants for 
the time evolution for < t < ti are displayed and Fig. 
|9]shows the coupling values for t\ < t < ti. Non-nearest- 
neighbor couplings are suppressed by 4 orders of mag- 
nitude due to the large distance of more than 260 /im 
between the corresponding ions. 



2. Second sequence: transforming a ID cluster state into a 
2D cluster state 

Now that a linear cluster state has been prepared dur- 
ing the first sequence of operations, qubits 1 and 4, 3 
and 6, and 5 and 8 need to be entangled in order to 
prepare a two-dimensional cluster state during the sec- 
ond sequence. This can be achieved in three steps. In 
the first step of the second sequence during time interval 
ti < t < t3, qubits 1 through 4 are stored in a common 
trap, e.g., with trap frequency u> = 2tt x 200kHz, while 
all other ions occupy single, individual traps with trap 
distances of 260 /im. Now, applying a magnetic-field gra- 
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4000 - 




FIG. 8: (Color online) Coupling constants for 8 ions stored in 
individual harmonic potential wells (shown in Fig. [7J. Each 
potential well is occupied by two ions, yielding uniform NN 
couplings of 3 kHz between ions number 1 and 2, 3 and 4, 
etc.. These coupling constants are required for the first time 
evolution in the first sequence of operations at the end of 
which a 2D cluster state is obtained. Non-nearest-neighbor 
couplings are suppressed by four orders of magnitude. 



4000^ 




FIG. 9: (Color online) After creating the NN couplings shown 
in Figure [8] the magnetic field gradient is switched off. Then, 
ions sharing a common oscillator are separated and trans- 
ported into four other potential wells, so that ions 2 and 3, 4 
and 5 etc. occupy now a common well. Switching on a mag- 
net field gradient of 100 T/m results in uniform NN couplings 
of 3 kHz again. This completes the second step of the first 
sequence of operations required for a 2D cluster. 



dient of 100 T/m results in a coupling of spins 1 and 4 of 
J14 = 1.24 kHz with time- evolution operator 

4 

U(tn) = exp(-- ^ 4,m*14^,»V)' ( 16 ) 
n,m— l;n<m 

Other couplings involving qubits 1 to 4 can be canceled 
by selective recoupling pulse sequences [34] . All couplings 
regarding qubits 5 to 8 can be neglected during this time 
interval, since those are approximately 0.5 Hz. The entire 
time evolution for coupling ions 1 and 4 including pulses 
to eliminate couplings involving qubits 2 and 3 is thus 
given by 

U t3 =U{t u /4)a X)2 U(t u /4)a x<2 x 

x (J x jU{t 14 / A)a x . 2 U {tu/ A)a x .2<J x ,z, 

where tu — an d o x ^ denotes the usual Pauli X op- 
erator acting on the Hilbert space of qubit k. To see that 
Ut 3 implements only the coupling between qubits 1 and 
4, it is convenient to think of the two rows as individual 
sequences. In each of the two sequences, all couplings 
involving qubit 2 are eliminated, because applying a x . 2 
before and after adds a minus sign to every a z . 2 

in this time-evolution due to the commutation relation 
of the Pauli matrices. Analogously, the a x ^ at the be- 
ginning and end of the second sequence affects nothing 
but adding a minus sign to all couplings involving qubit 
3. In total, all couplings to qubits 2 and 3 differ in their 
signs for exactly half of the time, and thus they are elim- 
inated. On the other hand, the coupling of qubit 1 to 4 



has the same sign during the entire sequence, such that 
U t3 = cxp(-if <r Zt ia zA ). 

At the end of time interval £14, the magnetic-field gra- 
dient is switched off, and the same procedure is repeated 
to entangle ions 3 and 6, respectively 5 and 8. Then 
the ions are arranged in a configuration, such that ions 3 
to 6 occupy a common potential well and all other ions 
are located in separate trap potentials. Switching on the 
magnetic-field gradient again results in 

i 6 

U(t 36 ) =exp(-- ^ Jn,mtze,Vz,nVz,m), (18) 
n,m— 3;n<m 

and the time evolution including recoupling pulses reads 
U u =U(t 36 /4)a X AU(t 36 /4)c7 xA x 

X (T x .$U{t 3 Q/A)a xA U{t 3(i /A)<T xA <T x £, 

with £ 36 = . Repeating this procedure for entangling 
ions 5 and 8, one obtains the following time evolution: 

U t5 =U(t 58 /4)a xfi U(t 5S /4)a xS x 

x (T x jU(t 5S /4:)a x . 6 U(t 58 /A)a Xt6 a x . 7 , 

where ins = 7^— and 

8 

U(t 58 ) = CXp(-- ^ Jn, m h 8 <7z,n<7z,m)- (21) 
n,m—5;n<m 
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The result of this procedure is the cluster state based on 
the graph shown in Fig. [2] 



I*) =u u u u u u u ta u tl \+) 9 

= exp(-i^(a Zyl a Z A + <J z . 3 a zfi - 

7 

+ ^z.^zfi + / "Vz,»<7z,j+i)h 



(22) 



This scheme can be scaled to any n x 2 cluster, which 
may then serve to simulate a n x m cluster. 



3. Summary of transport scheme for generating 2D cluster 

states 

A summary of the scheme that makes use of transport 
of ions over small distances is provided in Table [V] If 
we restrict the scheme to adiabatic ion transport, the 
required time tx for transport must obey tx » 2tt/v\. 
For v\ « 2tt x 200 kHz, we estimate tx = 50 /is which is 
still more than 1 order of magnitude less than the time 
required for entangling the qubits. The ions need only 
be linearly transported over distances of the order of the 
size of two electrode segments; in our concrete example 
considered above this amounts to a distance of around 
200 /im. 

Adiabatic switching of the magnetic-field gradient 
would require a time scale of order ts >> ^7, if only 
the qubit states |0) and |1) were present. However, usu- 
ally other nearby ionic states have to be included in these 
considerations. For example, in the case of 171 Yb + ions 
and choosing |0) = \Sy 2 ,F = 0) and |1) = \S 1/2 ,F = 
l,mp = 1) one would require ts >> ^ to avoid tran- 
sitions between Zeeman states, where A indicates the 
Zeeman splitting of the (S1/2, F = 1) state. In order 
to avoid zero crossings of the states, while the magnetic- 
field gradient is changed, a constant offset field is used, 
e.g., such that A = 27rlOMHz without gradient, so that 
we estimate ts ~ 1 A*s. 

During this switching process the qubits' phases are 
affected due to the spin-spin coupling and the evolution 
of the Zeeman state. Applying the decoupling scheme 
described in |34j with appropriate time intervals removes 
changes in the qubits' phases due to this Zeeman evolu- 
tion together with the undesired spin-spin couplings as 
described above. 



III. CONCLUSION 

We propose methods to generate cluster states of 
trapped ions confined in state-of-the-art segmented lin- 
ear ion traps by engineering their spin-spin coupling con- 
stants. Based on the idea of simulating a n x m cluster by 
a n x 2 cluster within the one-way model of quantum com- 
puting, we examined in Sec. II A a method previously 



suggested |T7j to prepare n x 2 clusters, and a novel idea 
based on creating individual potential wells for each ion. 
In addition, the superposition of harmonic potentials is 
discussed in order to engineer J couplings fulfilling suit- 
able periodicity relations to create small cluster states in 
one time-evolution step. In order to achieve sizeable cou- 
pling constants (in the kilohertz range) , all these methods 
require control over local electrostatic potentials with a 
spatial resolution of the order 10 /im - a typical interion 
distance in usual Paul traps. Therefore, it is of interest 
to investigate how suitable a trap with larger electrode 
structures is for generating 2D cluster states. 



In Sec. |IIB| a scheme for preparing n x 2 clusters is de- 
scribed for which larger trap electrode structures (of the 
order 100 fim) are sufficient. Here, nx2 clusters are pre- 
pared by first creating a linear cluster and subsequently 
enabling third- neighbor couplings. We showed that the 
generation of the linear cluster state can be accomplished 
with modern segmented ion traps by locating pairs of 
two ions in common harmonic oscillators, thus result- 
ing in uniform NN couplings. After separating the ions 
and subsequently merging them with the other nearest- 
neighbors, the second half of NN couplings are realized. 
The required third-neighbor couplings are achieved via 
selective recoupling techniques. 

Common to all schemes described in this paper is that 
entanglement is achieved solely by controlling dc voltages 
and currents: no coherent interaction between laser light 
and trapped ions is required. Refocusing pulses applied 
to individual ions consist of radio-frequency or microwave 
radiation depending on the choice of qubit [2"2l |2"51 |4"21 143] . 



The method described in Sec. |IIB| is explicitly worked 
out for the generation of a eight-qubit 2D cluster state. 
This scheme is also applicable to the generation of n x 2 
cluster states where n = k x 4, k = 2, 3,4, .... The recipe 
would then be to create in parallel k 2D cluster states 
of size 4x2. Then in, one additional step, ions at the 
edges of neighboring 4x2 clusters are entangled. For 
example, generating a 8 x 2 cluster proceeds as follows: 
first create two 2D cluster states of size 4x2, then by 
combining ions number 7, 8, 9, and 10 in one potential 
well, simultaneously entangle ion pairs 8/9 and 7/10 to 
complete a 16-qubit 2D cluster state. 

The physical arrangement of a 2D cluster state of size 
n x 2 is not required to be a linear ion string as was 
assumed so far. For instance, eight ions could reside in 
one area of a large 2D trap array |44j and communication 
between different areas, and thus entanglement, could be 
achieved by shuttling only the ions at the ends of each 
ion string. 

The spin-spin coupling that is used here does not re- 
quire cooling of the ion string to its motional ground 
state. Detailed calculations show that cooling to the 
Dopplcr limit is sufficient for suppressing unwanted ef- 
fects of thermal motion. This is beyond the scope of this 
paper and will be subject of a separate publication. 
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Step 


Operation 


1 


Transporting ions 1/2, 3/4, 5/6, 7/8 in common trap potentials 




Entangling ions 1/2, 3/4, 5/6, 7/8 


2 


Transporting ions 2/3, 4/5, 6/7 in common trap potentials 




Entangling of ions 2/3, 4/5, 6/7 


3 


Recombination of ions 1-4 and 5-8 




Entangling ions 1 and 4 using selective recoupling 


4 


Recombination of ions 3-6 




Entangling ions 3 and 6 using selective recoupling 


5 


Recombination of ions 5-8 




Entangling ions 5 and 8 via selective recoupling 



TABLE V: Summary of transport scheme for preparing a n x 2 cluster state in a segmented ion trap with magnetic-field gradient 
of 100 T/m. Steps 1 through 2 serve to create a linear cluster state, while steps 3 through 5 create additional third neighbor 
couplings that turn the one-dimensional cluster into a two-dimensional cluster. Here, the coupling constants Ji,j+i = 3.0 kHz 
for preparation of the linear cluster state which yields the gate times tj — 2 j^- +1 ~ 0.52ms, j — 1,2. In the case of 4 ions 
sharing a common potential well Ji,4, ^3,6, ^5,8 are given by 1.2kHz and we have tj = 2 j- 7r +1 = 1-3 ms, j = 3,4,5. The time 
scale for adiabatically transporting ions is tr 3> 2n/ui « 5 /is. The time scale for adiabatically turning on and off the magnetic 
field is tB 3> 2tt/A. When using m Yb + ions the relevant A indicates the Zeeman splitting of the (5*1/2, F — 1) state and 2tt/A 
is typically of order 0.1 /its 
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